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Abstract 

A thermodynamic formulation for moving granular material is proposed. The fluctuations 
due to the constant flux and dissipation of energy are controlled in a "granular" ensemble 
by a pressure p ("compression") which is conjugate to a contact volume ("contactopy"). 
The corresponding response function ( "dissipativity" ) describes how dissipation increases 
with p and should serve to identify the fluidization transition and 1/f noise. In the granular 
ensemble one can consider the granular medium as a gas of elastically colliding particles 
and define a "granular" temperature and other standard thermodynamic quantities. 
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Granular materials, like sand or powders, subjected to an external force will lo- 
cally perform rather statistical motion due to the random nature of the size and shape 
of grains and their contacts. One example is the motion of sand on a vibrating plate, say 
a loudspeaker [1-6]. At sufficiently high frequency the individual grains chaotically jump 
up and down forming a gas-like cloud of colliding particles. Other examples are displace- 
ments inside a shear-cell[7-ll] or flow down an inclined chute[ll-15] where in addition to 
a laminar flow with a well defined (average) velocity profile one has Brownian-like motion 
of the particles perpendicular to the flow direction. 

The above observations have inspired several authors to use ther mo dynamical con- 
cepts to describe granular media. On one hand a "granular temperature" T gr has been 
defined[7,16,17] as T gr = (v 2 ) — (v) 2 : i.e. proportional to the kinetic energy surplus with 
respect to the global motion. This temperature has been determined numerically as a 
function of various external parameters and material constants and under certain condi- 
tions consistency with experimental measurements was confirmed [7]. The drawback of the 
above definition is that it is only ther mo dynamically justified if an equipartition theorem 
exists which is not the case for granular particles since they dissipate energy at collisions. 

Edwards and collaborators [18-20] have put forward an entirely different, original idea: 
Based on the important observation that granular materials do not conserve energy while 
the entropy S is well defined they proposed to consider the volume V to replace the 
internal energy in the usual thermodynamic formalism. In this way a temperature-like 
quantity X = dV/dS which they called "compactivity" can be defined. Although formally 
correct, this formalism is not easy to justify from first principles. In particular, in many 
real situations like on the vibrating table or on an inclined plane, the volume is not well 
limited at large heights. While Edwards's approach seems intuitively correct for dense 
packings and the definition of T gr reasonable in the limit of strong internal motions or 
weak dissipation they fail in the corresponding opposite limit. 

Purpose of the present note is to propose a different thermodynamic approach to 
granular materials founded on similar principles as equilibrium thermodynamics and which 
should at least partially incorporate the intuitive pictures of previous work. 

As opposed to usual thermodynamics of molecular gases the elementary units of gran- 
ular materials are mesoscopic grains consisting of many atoms each (10 15 — 10 25 ). When 
these object interact (collide) the Lennard- Jones potentials of the individual atoms are 
unimportant and completely different mechanisms must be considered. It is important 
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that on a microscopic scale the surface of the grains is rough. Solid friction is the immedi- 
ate consequence: When two touching grains are at rest with respect to each other a finite 
force F s is needed to trigger relative motion (static friction), while moving against each 
other a finite force Fa is needed to maintain the motion (dynamic friction). Fa < F s and 
both only depend on the normal force and neither on the velocity nor on the area of contact 
(Coulomb law). No doubt, this picture is idealized and an entire discipline, called tribol- 
ogy, has evolved to study solid friction in depth[21]. For our purpose it is, however, more 
convenient to remain on the simple text-book level. The solid friction has the crucial con- 
sequence that on the level of the elementary units, namely the grains, the system does not 
conserve energy as opposed to molecular thermodynamics. Another source of dissipation 
can be plastic deformation of grains due to the normal force acting at collisions. 

If energy is not constantly pumped into a granular system it will stop moving and fall 
into one of its static configurations. Constant motion of the grains can only be produced 
when there is a steady state of energy flux. We are, however, not interested in this flux 
itself also because it is difficult to measure experimentally. We just want to describe the 
motion of the granular particles in a similar way as one describes the motion of molecules 
in a gas at a given temperature. The presence of the energy flux and the fact that on the 
level of the grains on which we want to formulate a thermodynamics the energy is locally 
dissipated (i.e. not conserved) will, however, force us to introduce concepts beyond that 
of usual equilibrium thermodynamics. 

We will assume typical conditions for local "equilibrium": Most experiments have 
velocity and density gradients [1-1 5] and in those cases only a subsystem spatially small 
compared to the gradient should be considered. An eventual energy flux into the system 
should distribute the energy over it homogeneously. This constraint can also reduce the size 
of the subsystem. Outside this subsystem a generalized "heat bath" is assumed. Spatial 
and temporal averages should be exchangeable ( "ergodicity" ) . We will in fact in the 
following consider temporal averaging for practical (numerical) purposes. The averaging 
procedure can even be complicated [4,6] by the existence of density waves. 

It is important to notice that the dissipated energy is of course only lost on the 
mesoscopic level - microscopically this energy will be transformed essentially into heat 
and blown away by the surrounding air. This gives us a reasonable starting point for the 
formulation of an analogy to usual thermodynamics. It seems natural to consider energy 
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conservation as the first "thermodynamical principle" : 

AI = AE int + AD . (f) 

The internal energy Ei nt is like in traditional thermodynamics the kinetic and potential 
energy of all the degrees of freedom of the grains as elastic bodies (translation, rotation, 
elasticity, etc). AD is the energy dissipated in a given time and AI is the energy that 
was pumped into the system while AD was dissipated in order to maintain the steady 
stated Usually AI is some kind of work (gravity on the inclined plane, \Auj 2 on the 
loudspeaker, etc). If one allows for changes in the volume of the system then eq. (1) will 
become AI — AE int + AD + AW where AW is the work done to change the volume. 

Let us give to the excess dissipated energy AV = AD — AI in the following the 
nick-name "dissipate" . We will deal with V in a similar way as one treats in usual ther- 
modynamics the heat. Like the heat, the dissipate is not a potential since it depends on 
the process by which a given state is reached. It does, however, not stem from the kinetic 
energy of the particles as the heat in a molecular gas but is due to collisions, i.e. two 
particles coming together, touching and separating again. 

The dissipated energy is proportional to the sum of normal forces f n that push the 
particles together during collision i. One can therefore express changes in T> as 

SV = p5C (2) 

where p is an internal pressure acting at collisions that we shall call "compression"^ It 
can be defined as p = p{f! n /A i ) where f n is the normal force and Ai the area of contact 
of collision i and the average is performed over all collisions, p is the density of collisions, 
defined as the number of collisions per unit volume and unit time. It is easy to determine 
p numerically. When the particles do not have collisions the compression is zero and no 
energy is dissipated. The quantity C, which we will in the following call "contactopy" , in 

t In contrast to traditional thermodynamics we have an energy flux and the dissipated 
energy itself increases with time. Therefore one can formulate eq. (1) alternatively as 
Ji = E int + D where Jj is the energy flux into the system. Although in some cases this 
description gives a more intuitive physical picture we will prefer in the following to argue 
only in terms of changes AD and AI during a fixed time interval. 

■f In real collisions also shear forces can contribute to dissipation so that p would then 
be a like a stress tensor. For that reason we explicitely did not call it pressure. 
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analogy to the entropy, has the dimension of a volume (contact volume). It is defined as 
the conjugate variable to the compression p. 

The contactopy has contributions due to plasticity and due to dynamic and static fric- 
tion. Let us in the following argue for a geometrical interpretation of C and consider first 
the two contributions from friction. The dynamic (or better kinetic) part of the contactopy 
is proportional to £\ A$i where ti is the distance over which two solid grains slip during 
collision i. Since li is given by the collision time multiplied by the velocity of the particles 
this part is proportional to the particle overlap volume V ov that one has (for technical 
reasons) in molecular dynamic simulations [4, 6, 10, 15, 2 2]. Apart from geometrical prefac- 
tors the proportionality constant is the dimensionless dynamic friction coefficient [id, i.e. 
a material constant of the grains. The overlap volume V ov can be defined more precisely 
as the sterically excluded volume that would arise if the centers of mass of the particles 
follow the real trajectories but one does not take into account the elasto-plastic deforma- 
tion. The static (or potential) contribution of friction to the contactopy C only comes 
into play when the elastic (potential) energy of two unlocking particles that were sticking 
is released. It depends on the penetration depth d{ at collision i because this determines 
the amount of material that will be compressed (or fragmented). Therefore this second 
contribution is proportional to ^ i-e. also proportional to V ov . The proportionality 
factor contains the static friction coefficient /j, s and the Young modulus Y of the grains. 
The other contribution to the potential part of the contactopy comes from plasticity and 
is proportional to the size of the plastic zone, i.e. again to V ov , when the material de- 
pendent plastic yield force F p is reached. The complex stick-slip mechanism[10,23] that is 
triggered between two rigid grains by the interplay of static and dynamic friction makes it 
difficult to determine precisely the resulting material dependent constant 7(^5 Us, F P ,Y) 
that following the simple arguments given above relates via 

C = ^V ov (3) 
P 

the contactopy to the (average) overlap volume V ov . The division by the number p' of 
collisions per unit volume in eq. (3) assures that C is an extensive quantity. The propor- 
tionality of eq. (3) can be checked numerically by measuring independently the input of 
energy and the average internal pressure to get C and summing up the overlap volumina 
of the collisions to get V ov . 

The contactopy plays a central role here and in similarity to the approach of Ed- 
wards [18] can be interpreted as a volume. But although it ressembles the internal energy 



5 



in having a kinetic and a potential contribution it is analogous in our thermodynamic 
formalism to the entropy. Because of eq. (3) the contactopy represents a geometric charac- 
terization of the system. This makes it likely to be a total differential dC, i.e. independent 
on how the system was driven into its state, while in contrast, the dissipate depends on 
the work done on the system. Numerically this could be checked by monitoring V ov but 
the ultimate test should be experimental in analogy to Carnot's experiments. 

The idea to define a potential for a dissipative system actually dates back to Lord 
Rayleigh[24] ("dissipation function") and has been worked out in detail by Prigogine and 
collaborators [2 5]. The contactopy is on one hand a concrete example for such a potential 
on the other hand it does not only contain the statistical aspects of an "internal entropy 
production" [25] but has for physical reasons the dimension of a volume. It would in fact 
be important to work out a statistical interpretation of C in the (space-time) phase space 
of the collision events in analogy to Boltzmann's statistical definition of the entropy in the 
space of all configurations of positions and velocities. 

The "equilibrium" - which is in fact a steady state driven by the energy flux - can now 
be defined as the ensemble minimizing at fixed E int the contactopy (instead of maximizing 
the entropy). One can postulate an analog to the second law of thermodynamics that any 
change of state at constant internal energy E int should decrease the contactopy 

AC < . (4) 

Physically such a behaviour seems naturally be driven by the elastic repulsion between col- 
liding (overlapping) grains and the tendency of the system to prefer many smaller collisions 
to a few strong ones. The third law of thermodynamics, namely that at zero compression 
there is no overlap between grains, i.e. vanishing contactopy, is less evident. Numerical 
tests of the above statements should be performed taking into account that as mentioned 
before they are valid in (sub) systems into which the energy flow allows for homogeneous 
dissipation. 

As in usual thermodynamics one can now work in different ensembles. One can fix 
either the compression p which we shall call the "granular ensemble" or the contactopy 
(let's call it "atomistic ensemble"). Since in practice (experimentally and numerically), 
however, the later case is difficult to implement we will in the following usually consider the 
granular ensemble. On top of this we can build up the traditional body of thermodynamics 
as if the grains were a gas of particles interacting elastically. We can fix or free the number 
N of particles, define a "granular" temperature T g and entropy S or impose to the system 
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either an external volume V or an external pressure p. A novelty for granular media is 
that one could also impose an external shear r or its "conjugate", the dilatancy Vd[26]. 
A granular potential G r can be defined as 



G r = E int + pC (5) 

which depends on p and the extensive variables N, V, S and Vd- An immediate consequence 
of eq. (4) is that at constant compression p the equilibrium is given by the minimum of 
G r . The analog to the specific heat might be called "dissipativity" k defined as 

dp dp 

This is a new quantity characterizing the granular medium which measures how much more 
energy can be dissipated if the compression is increased. It could be measured directly by 
numerically evaluating the derivative of eq. (6) or through the fluctuations of the energy 
in a thermally closed system, i.e. surrounded by (infinitely) heavy and stiff walls. The 
dissipativity k should be positive and go to zero for p — > and p — > oo. Interesting 
for practical purposes is that k contains through the 7 of eq. (3) the material dependent 
properties concerning friction, among others also the stick-slip mechanism between grains. 
If a fluctuation-dissipation theorem for the response function k is valid then one might 
identify 1/ / noise[23,27] from its frequency dependence over time scales proportional to 
the size of the grains or even over larger time scales when collective phenomena like arching 
or bridge-collapsing[8,9] come into play. 

We know that there exists a "fluidization transition" in granular media between a 
regime of block motion at low energy flow to a gas-like collisional regime at high energy 
flow [3, 4]. This transition could be driven by changing p: For small p the potential part 
of C dominates (block motion) and for large p the kinetic part of C is relevant (colli- 
sional motion). It seems likely that the transition point is given by a singularity of the 
dissipativity k. This could be checked experimentally and numerically. 

In the granular ensemble a "dissipate" bath (instead of a heat bath) is coupled to the 
system and consequently the internal energy E int of the granular material is a fluctuating 
quantity. In order to give Boltzmann's statistical interpretation to the entropy S it is 
therefore conceptually better to work in the atomistic ensemble: A "state" is given by 
the positions, orientations, linear and angular velocities of the grains as rigid bodies. In 
fact the entropy is well-defined as noted already in ref. 18. A reasonable definition for a 
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granular" temperature T g would then be: 




( 



dG r 



dS 



) 



p 



(7) 



Note that it is similar but not identical to the granular temperature defined by previous 
authors[7,16,17]. T g is the variable that controls the granular canonical ensemble with the 
granular free energy F g as potential, defined as the Legendre transformation of the granular 
potential: F g = G r —T g S . In equilibrium F g should have a minimum. The usual Boltzmann 
distribution determines the statistical weights of the states in this ensemble. As in the case 
of the usual temperature one can measure T g by monitoring the exchange of internal energy 
between a subsystem and its heat bath which should obey this Boltzmann distribution. 
Also a direct measurement of T g by changing an external pressure (see eg. (17.1) of ref. 28) 
should be possible. A specific heat can be defined as a derivation with respect to T g . 

Experimentally p and T g are independent control parameters of the system: Since T g 
increases with the kinetic energy of the particles it is essentially controlled by the amount 
AI of energy that is fed into the system per unit time. The compression p or better 
the quantity tJ?- also depends on the density of collisions and can therefore increase by 
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fragmenting the grains into smaller pieces. (Note that when a given grain is split into 
eight pieces, the cross section of each individual piece decreases by a factor four, so that p 
will increase by two.) One can therefore, by changing the grain diameters and AI modify 
-^r- and T g independently and can therefore consider phase diagrams in the ^- - T g plane. 

As already mentioned one can also go to other ensembles by changing variables via 
further Legendre transformations. One can liberate the number of grains and introduce 
a granular grandcanonical potential controlled by a chemical potential. One can also fix 
an external pressure and calculate the average volume. The fact that both volume and 
contactopy will then be conjugate to the pressure could explain why in granular media 
one finds in equilibrium macroscopic density fluctuations [29] as opposed to usual fluids 
or gases. Of practical interest is also to fix an external shear and measure the average 
expansion, i.e. the dilatancy[8,9,26]. 

It is useful to note that in the case when friction and plasticity vanish the system 
does not dissipate energy anymore, the contactopy will be zero and G r = Ei nt . In that 
case the atomistic and granular ensembles are identical and classical thermodynamics is 
immediately recovered. Our formalism is therefore a genuine generalization of equilibrium 
thermodynamics. 
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We have described within a thermodynamic formalism the fluctuations arising from 
the constant flux and dissipation of energy that drives a granular material's kinematic 
behaviour. By separating the dissipative degree's of freedom (friction and plasticity) from 
the conservative ones (translation, rotation, elasticity) we define a "granular ensemble" 
coupled to a "dissipate bath" which is in fact the one in which experimental and numerical 
measurements are usually performed. We introduce a potential that we call "contactopy" 
and argue that it is proportional to the steric overlap volume of the collisions which the 
particles would have had per time unit if while following the real trajectories they had no 
elasto-plastic deformation. It would be interesting to give also a statistical interpretation to 
the contactopy in order to define it as a dissipative potential [2 5]. The fluctuating internal 
energy is replaced by a granular potential controlled by an intensive variable that we call 
"compression", which is conjugate to the contactopy. Going into a granular canonical en- 
semble we define a "granular temperature" similar to the one defined previously [16, 17]. We 
propose various numerical and experimental tests for the assumptions that we have made 
in our theory and suggest that a frequency dependent "dissipativity" should characterize 
the stick-slip behaviour of the material and the transition to fluidization. 

I thank R. Balian, J. Kertesz, M. Kiessling, J. Lebowitz, C. Moukarzel, S. Roux, D. 
Stauffer, D. Wolf, R.K.P. Zia and in particular H.G. Herrmann for encouraging discussions. 
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